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The proton spin puzzle has challenged our understanding of QCD for the last 20 years. 
We survey new developments in theory and experiment. The proton spin puzzle seems to 
C^h' be telling us about the interplay of valence quarks with chiral dynamics and the complex 

' vacuum structure of QCD. 
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1. Introduction 



Protons behave like spinning tops. Unlike classical tops, however, the spin of these 
particles is an intrinsic quantum mechanical phenomenon. This spin is responsible 
for many fundamental properties of matter, including the proton's magnetic mo- 
ment, the different phases of matter in low-temperature physics, the properties of 
neutron stars, and the stability of the known universe. How is the proton's spin 
built up from the spin and orbital angular momentum of the quarks and gluons 
inside ? 

Polarized deep inelastic scattering experiments have revealed a small value for 
the nucleon's flavour-singlet axial-charge g^| P Dis ~ 0.3 suggesting that the quarks' 
' intrinsic spin contributes little of the proton's spin. The challenge to understand the 

spin structure of the protorP^ has inspired a vast programme of theoretical activity 
and new experiments. Why is the quark spin content g^lpDis so small ? 

We start by recalling the g\ spin sum-rules, which are derived starting from the 
dispersion relation for polarized photon-nucleon scattering and, for deep inelastic 
scattering, the light-cone operator product expansion. One finds that the first mo- 
ment of the gi spin structure function is related to the scale-invariant axial-charges 
of the target nucleon by 
1 

dx g^(x,Q 2 ) = 

o 



(1) 
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Here , g^ an d fl^i linv are the isovector, SU(3) octet and scale-invariant flavour- 
singlet axial-charges respectively. The flavour non-singlet c^se and singlet cse Wilson 
coefficients are calculable in ^-loop perturbative QCDP The term /3oo represents 
a possible leading-twist subtraction constant from the circle at infinity when one 
closes the contour in the complex plane in the dispersion relation^ If finite, the 
subtraction constant affects just the first moment. The first moment of 171 plus the 
subtraction constant, if finite, is equal to the axial-charge contribution. 
In terms of the flavour dependent axial-charges 

2Ms, l Aq = (p, s| 57^75 <?|p, s) (2) 

the isovector, octet and singlet axial-charges are: 

g^ = Au-Ad 

g^ = Am + Ad - 2As 

gf\ im /E(a s ) = gf = Au + Ad + As. (3) 

Here E(a s ) = exp f™ s da s 7(<S s )//3(a ; s ) is a renormalization group factor which cor- 
rects for the (two loop) non-zero anomalous dimension 7(a s ) of the singlet axial- 
vector current P = ^7^751* + dj^j^d + sj^j^s , and which goes to one in the 
limit Q 2 — > 00; (3(a s ) is the QCD beta function. The singlet axial-charge, g^linv, is 
independent of the renormalization scale (i and corresponds to g^ (Q 2 ) evaluated 
in the limit Q 2 — > 00. The axial-charges g^ and g^ are renormalization group 
invariants. 

If one assumes no twist-two subtraction constant (/Jco = 0(1/Q 2 )) the axial- 
charge contributions saturate the first moment at leading twist. The isovector axial- 
charge is measured independently in neutron /3-decays {g^ = 1.270 ± O.OOcP) and 
the octet axial-charge is commonly taken to be the value extracted from hyperon 
/3-decays assuming a 2-parameter SU(3) fit (g^ — 0-58 ± 0.03^). Using the sum- 
rule for the first moment of gi, given in Eq. (1), polarized deep inelastic scattering 
experiments have been interpreted in terms of a small value for the flavour-singlet 
axial-charge. If we take g^ = 0.58 ±0.03, then inclusive g\ data with Q 2 > 1 GeV 2 
give^ 

9 { a l P Dis,Q»-foo = 0.33 ± 0.03(stat.) ± 0.05(syst.) (4) 

- considerably smaller than the value of g A quoted above. 

In the naive parton model g^| P Dis is interpreted as the fraction of the proton's 
spin which is carried by the intrinsic spin of its quark and antiquark constituents. 
When combined with g^ = 0.58 ±0.03 the value of <?^| p dis in Eq.(4) corresponds 
to a negative strange-quark polarization 

As Q 2^ oc = i(gi 0) | P Dis,Q^co - ffi 8) ) = -0.08 ± 0.01 (stat.) ± 0.02(syst.) (5) 

- that is, polarized in the opposite direction to the spin of the proton. 

What physics separates the values of the octet and singlet axial-charges ? 
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2. Spin and the singlet axial-charge g^ 1 

There are two key issues: the physics interpretation of the flavour-singlet axial- 
charge and possible SU(3) breaking in the extraction of from hyperon 
/3-decays. 

First consider g^\ Gluonic information feeds into g^ through the QCD axial 
anomaly. QCD theoretical analysis leads to the formulsP^Hni 

^=(£A<Z-3|W) +C00. (6) 

V q / partons 

Here Ag partons is the amount of spin carried by polarized gluons in the polarized 
proton (a s Ag ~ constant as Q 2 —> ocP^) anc [ Aq pm - tons measures the spin carried by 
quarks and antiquarks carrying "soft" transverse momentum k 2 ~ P 2 , m 2 where P 
is a typical gluon virtuality and m is the light quark mass. The polarized gluon term 
is associated with events in polarized deep inelastic scattering where the hard photon 
strikes a quark or antiquark generated from photon-gluon fusion and carrying k 2 ~ 
Q 2 ! 10 * 11 ! Coo denotes a potential non-perturbative gluon topological contribution 
which is associated with the possible subtraction constant in the dispersion relation 
for g x and Bjorken x = OP g { ^> | pD is = 9a ~ C °°- 

The subtraction constant, if finite, is a non-perturbative effect and vanishes in 
perturbative QCD. It is sensitive to the mechanism of axial U(l) symmetry break- 
ing and the realization of axial U(l) symmetry breaking by instantons: spontaneous 
U(l) symmetry breaking by instantons naturally generates a subtraction constant 
whereas explicit symmetry breaking does notP^l The QCD vacuum is a Bloch su- 
perposition of states characterised by non-vanishing topological winding number 
and non-trivial chiral properties. When we put a valence quark into this vacuum it 
can act as a source which polarizes the QCD vacuum with net result that the spin 
"dissolves" and some fraction of the spin of the constituent quark is associated with 
non-local gluon topology with support only at Bjorken x = 0. 

Possible explanations for the small value of g^lpDis extracted from the polar- 
ized deep inelastic experiments include screening from positive gluon polarization, 
negative strangeness polarization in the nucleon, a subtraction at infinity in the 
dispersion relation for g\ associated with non-perturbative gluon topology and con- 
nections to axial U(l) dynamics JA3HHD as we \\ as possible SU(3) breaking in g^ - 
possibly as large as 20%P^^ The QCD axial anomaly decouples from the non- 
singlets and ■ 

One would like to understand the dynamics which appears to suppress the singlet 
axial-charge extracted from polarized deep inelastic scattering relative to the OZI 
prediction g^ = g^ and also the sum-rule for the longitudinal spin structure of 
the nucleon 

\ = \Y.^ + ^9 + L q + L g (7) 
where L q and L g denote the orbital angular momentum contributions. There is 
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presently a vigorous programme to disentangle the different contributions. Key ex- 
periments include semi-inclusive polarized deep inelastic scattering (COMPASS and 
HERMES) and polarized proton-proton collisions (PHENIX and STAR at RHIC), 
as well as deeply virtual Compton scattering to learn about total angular momen- 
tum. 

3. The shape of g\ 

To understand the proton spin puzzle, it is interesting to look at the x dependence 
of the measured g\ spin structure function. Deep inelastic measurements of g\ have 
been performed in experiments at CERN, DESY, JLab and SLAC. There is a general 
consistency among all data sets. COMPASS are yielding precise new data at small 
x, down to x ~ 0.004. JLab are focussed on the large x region. 

Precise measurements of the deuteron spin structure function gf show the re- 
markable feature that gf is consistent with zero in the small x region between 0.004 
and 0.02P In contrast, the isovector part of g\ is observed to rise at small x as 
~ 2;-0-22±o.o7 anc j j g mucn bigger than the isoscalar part of #iP^ This is in sharp 
contrast to the situation in the unpolarizcd structure function F% where the small x 
region is dominated by isoscalar pomeron exchange. The g\~ n data are consistent 

with quark model and perturbative QCD predictions in the valence region x > 0.2PS1 

(3) 

The size of g A forces us to accept a large contribution from small x and the ob- 
served rise in g\~ n is in excellent agreement with the prediction gf n ~ a; -0 22 of 
hard Regge exchange - in particular a possible a\ hard-pomeron cut!^ involving the 
hard-pomeron which seems to play an important role in unpolarized deep inelastic 
scattering P2l 

The "missing spin" is associated with a "collapse" in the isosinglet part of g\ 
to something close to zero instead of a valence-like rise for x less than about 0.02. 
This isosinglet part is the sum of SU(3)-flavour singlet and octet contributions. If 
there were a large positive polarized gluon contribution to the proton's spin, this 
would act to drive the small x part of the singlet part of g\ negative 2 ^ - that is, 
acting in the opposite direction to any valence-like rise at small x. However, gluon 
polarization measurements at COMPASS, HERMES and RHIC constrain this spin 
contribution to be small in measured kinematics meaning that the sum of valence 
and sea quark contributions is suppressed at small x. (Soft Regge theory predicts 
that the singlet term should behave as ~ TVlnx in the small x limit, with the 
coefficient N to be determined from experiment ! 24 l 25 l) 

There is presently a vigorous programme to disentangle the different contri- 
butions involving experiments in semi-inclusive polarized deep inelastic scattering 
and polarized proton-proton collisions pSHH3 These direct measurements show no 
evidence for negative polarized strangeness in the region x > 0.006 (in apparent 
contrast to the extraction of negative strangeness polarization extracted from in- 
clusive measurements of g\). For gluon polarization, present measurements suggest 
I - 3fjA5| < 0.06 corresponding to |A#| < 0.4 with a s ~ 0.3. That is, they 
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are not able to account for the difference (<7^| p dis — 9a^) ~ —0.25 obtained via 
Eq.(4). An independent measurement of the strange-quark axial-charge could be 
made through neutrino-proton elastic scattering P3 The axial-charge measured in 
vp elastic scattering is independent of any assumptions about the presence or ab- 
sence of a subtraction at infinity in the dispersion relation for g\ and the x ~ 
behaviour of g%. Further measurements to push the small x frontier in polarized 
deep inelastic scattering would be possible with a polarized ep collider.^ 

4. SU(3) breaking and 

Given that the contributions to g A from the measured distribution As and from 
— 3§^A<7 are small, it is worthwhile to ask about the value of g A . The value 0.58 is 
extracted from a 2 parameter fit to hyperon /3-decays in terms of the SU(3) constants 
F = 0.46 and D = 0.8(P- see Table 1. The fit is good to - 20% accuracyP^ The 

(8) 

uncertainty quoted for g A has been a matter of some debate. There is considerable 
evidence that SU(3) symmetry may be badly broken and some have suggested that 
the error on g^ should be as large as 25% More sophisticated fits will also include 
chiral corrections. Calculations of non-singlet axial-charges in relativistic constituent 
quark models are sensitive to the confinement potential, effective colour-hypcrfinc 
interactionpSEH pion and kaon clouds plus additional wavefunction corrections^ 
chosen to reproduce the physical value of g A . 

This physics has recently been investigated by Bass and Thomas^ within the 
Cloudy Bag model (CBM^SHH which has the attractive feature that when pion cloud 
and quark mass effects are turned off the model reproduces the SU(3) analysis. One 
finds that chiral corrections significantly reduce the value of g A . This, in turn, has 
the effect of increasing the value of g^ A | p dis an d consequently reducing the absolute 
value of the "polarized strangeness" extracted from inclusive polarized deep inelastic 
scattering. 

The Cloudy Bag^ was designed to model confinement and spontaneous chiral 
symmetry breaking, taking into account pion physics and the manifest breakdown 
of chiral symmetry at the bag surface in the MIT bag. If we wish to describe proton 
spin data including matrix elements of J^ 5 , J* 5 and J M 5, then we would like to 
know that the model versions of these currents satisfy the relevant Ward identities. 

(3) (8) 

For the scale- invariant non-singlet axial-charges g A and g A , corresponding to the 

Table 1. QaIqv from /3-decays with F = 0.46 and D = 0.80, together with 
the mathematical form predicted in the MIT Bag with effective colour-hyperfine 
interaction (see text and Ref. [32]). 



Process measurement SU(3) combination Fit value MIT + OGE 



n — 


' P 


1.270 ±0.003 


F + 


D 


1.26 


§B' 


+ G 


A 


-+p 


0.718 ±0.015 


F + 


& 


0.73 




B' 


£- 


— > n 


-0.340 ±0.017 


F - 


D 


-0.34 


l B ' 


- 2G 




->A° 


0.25 ±0.05 


F - 




0.19 


-G 


3° - 




1.21 ±0.05 


F + 


D 


1.26 


\B> 


+ G 



April 29, 2010 0:10 



WSPC - Proceedings Trim Size: 9.75in x 6.5in gm 



6 

matrix elements of partially conserved currents, the model is well designed to make 
a solid prediction. 

The effective colour-hyperfine interaction has the quantum numbers of one-gluon 
exchange (OGE). In models of hadron spectroscopy this interaction plays an im- 
portant role in the nucleon-A and E — A mass differences, as well as the nucleon 
magnetic moments^ and the spin and flavor dependence of parton distribution func- 
tions.^ It shifts total angular-momentum between spin and orbital contributions 
and, therefore, also contributes to model calculations of the octet axial-chargesf^^H 
In Bag model calculations one also needs to include wavefunction corrections asso- 
ciated with the well known issue that, for the MIT and Cloudy Bag models, the 
nucleon wavefunction is not translationally invariant and the centre of mass is not 
fixed. To compare the model results with experiment we take the view^l that, in 

(3) 

principle, the model - with corrections - should give the experimental value of g A . 
We therefore choose the centre-of-mass factor phcnomcnologically to give the ex- 

(3) 

pcrimental value of g A . This then fixes the parameters of the model and allows us 

(81 

to use it to make a model prediction for g A . 

Without pion cloud corrections the MIT Bag with centre of mass corrections 
reproduces the SU(3) analysis of the axial-charges extracted from /3-decays. This is 
illustrated in Table 1. Without additional physics input, e.g. pion chiral corrections, 
there is a simple algebraic relation between the SU(3) parameters F and D, the bag 
parameter B' and the OGE correction G: F = \B' - \G and D = B' + §G. The 
numerical agreement is very goodP^ 

The pion cloud of the nucleon also renormalizes the nucleon's axial-charges by 
shifting intrinsic spin into orbital angular momentumPS^Hln the Cloudy Bag Model 
( CBM) the nucleon wavefunction is written as a Fock expansion in terms of a bare 
MIT nucleon, |N), and baryon-pion, |Ntt) and |A7r), Fock states. The probabilities 
to find the nucleon in each Fock component are determined phenomenologically by 
fitting to a wealth of nucleon observablesP^ The expansion converges rapidly and 
we may safely truncate the Fock expansion at the one pion level. When we calculate 

(8) 

the pion and kaon cloud chiral corrections to g A we also have to choose the chiral 
representation, in particular whether to use the original surface coupling or the later 
volume coupling version of the Cloudy Bag model. 

(8) 

The extent of the reduction in g A depends upon the version of the CBM used, 
lying in the range 0.49 ± 0.02 for the original CBM and 0.42 ± 0.02 for the volume 
coupling version.^ These changes alone raise the value of g^ A | p dis,q 2 ^oo derived 
from the experimental data from 0.33±0.03(stat.)±0.05(syst.) to 0.35±0.03(stat.)± 
0.05(syst.) and 0.37 ± 0.03(stat.) ± 0.05(syst.), respectively. Both of these values 
have the effect of reducing the level of OZI violation associated with the difference 
.gi 0) | P Dis-si 8) from -0.25±0.07tojust -0.14±0.06 and -0.05±0.06, respectively. 
It is this OZI violation which eventually needs to be explained in terms of singlet 
degrees of freedom: effects associated with polarized glue and/or a topological effect 
associated with x = 0. 

The uncertainty in this model calculation lies in the small ambiguity between the 
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two chiral representations that one can choose. In order to quote an overall value 
that properly encompasses these possibilities we follow the Particle Data Group 
procedure^ for combining data that may not be compatible to estimate the overall 
error, finding a combined value of g A = 0.46 ± 0.05 (with the corresponding semi- 
classical singlet axial-charge or spin fraction being 0.42 ± 0.07 before inclusion of 
gluonic effects).^ Note that the error ±0.02 on g^ quoted for each version of the 
model follows from varying over the phenomcnological range of possible pion pa- 
rameters within first order perturbation theory. In terms of analogy to experimental 
errors, the ±0.02 is like a statistical error and the final ±0.05 error includes sys- 

(8"t 

tematic effects. With this final value for g A the corresponding experimental value 
of g^'lpDis would increase to g A | p dis = 0-36 ± 0.03 ± 0.05. 

5. Towards possible understanding 

Where are we in our understanding of the spin structure of the proton and the small 
value of g^ A | p dis ? Measurements of valence, gluon and sea polarization suggest that 
the polarized glue term — 3^-A<7 par t ons and strange quark contribution As par tons in 
Eq.(6) are unable to resolve the small value of g^| P Dis- Two explanations are 
suggested within the theoretical and experimental uncertainties depending upon 
the magnitude of SU(3) breaking in the nucleon and hyperon axial-charges. One 
is a value of g^ A ~ 0.5 (as suggested by the surface coupling model) plus an axial 
U(l) topological effect at x = associated with a finite subtraction constant in 
the g± dispersion relation. The second is a much larger pion cloud reduction of 
g A to a value ~ 0.4 (as suggested by the volume coupling model in first order pion 
cloud perturbation theory) . Combining the theoretical error on the pion cloud chiral 
corrections embraces both possibilities. The proton spin puzzle seems to be telling 
us about the interplay of valence quarks with chiral dynamics and the complex 
vacuum structure of QCD. 
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